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Abstract. In this paper we describe a renormalon- inspired model of twist-4 power corrections 
to the light-cone pion distribution amplitude, and compare it to the results obtained using the 
conformal wave expansion. We hnd that the overall functional form predicted in this renormalon 
model is similar to the one predicted in the standard approach to higher twist contributions 
involving an expansion in conformal operators. However, small discrepancies at the end-points 
allow for a discussion on the convergence properties of the conformal expansion. 



NBI-HE-99-33 



* andersen@nbi.dk 



1 Introduction 



Hadron distribution amplitudes play a key role in the study of hard exclusive processes in QCD. 
At large Q^, amplitudes for such processes can be written as a convolution of a process dependent 
hard scattering amplitude which is calculable in perturbation theory, and a process independent 
non-perturbative distribution amplitude, one for each hadron participating in the process 
The leading twist-2 distribution amplitude for the pion at asymptotically large is given by 

(piu) = 6u{l-u), (1) 

to which there are corrections at finite in the form of Gegenbauer polynomials. It is, however, 
generally believed that there will be large power suppressed corrections to the usual factorisation 
at accessible energies. For some observables (eg. the pion transition form factor for the process 
vr 7*7* and light-cone sum rules for the decay B — > vrei^), these power suppressed contributions 
are calculable in terms of process independent higher twist distribution amplitudes coming from 
higher Fock states and transverse momentum of the quarks. However, for eg. the pion form 
factor one has to also take into account so-called soft contributions (see eg. |§, ^) that generally 
result in power corrections to the process dependent hard scattering amplitudes. This paper will 
not consider the problem of factorisation at the power suppressed levels or the soft contributions. 
Instead, we will focus on the higher twist distribution amplitudes from transverse momentum 
and Fock states with different spin components. 

A gauge-invariant way to study the expansion of the pion distribution amplitude in terms of 
increasing twist is to consider the vacuum-to-pion matrix element. We parametrise it as Q 
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where (t>{u) is the asymptotic leading twist contribution (|l|) and [0, x] denotes the path-ordered 
gauge factor along the straight line between and x 



[x,y] = Pexp 



ig / dt{x - y)f,A'' {tx + (1 - t)y) 







(3) 



gi{u) and 52 (^) parametrise the twist-4 part and correspond to the contributions from the 
transverse momentum of the quarks and from states of the quarks with a different spin projection 
'ip'y-i/j, respectively. 

The standard procedure for estimating higher twist contributions is to exploit a connection 
between gi{u), g2{u) and the three-particle amplitudes. This connection is derived using the 
QCD equations of motion and allow one to calculate gi{u) and 52 (^) through an expansion in 
conformal operators of the three-particle operator. It is generally believed that the anomalous 
dimensions of operators increase with increasing conformal spin, and so only a few terms in this 
expansion should be needed to get reliable results. The analysis in Q includes the contribution to 
twist-4 distribution amplitudes from the operators of lowest and next-to lowest conformal spin. 
In this framework there has been recent progress in the study of pre- asymptotic behaviour of 
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the meson distribution amplitudes, extending the analysis to vector mesons and including mass- 
corrections p. The results are being extensively used, e.g. in light-cone sum rule calculations 
of B-meson decay form factors etc. 

In the following we will describe a different approach that does not rely on the hierarchy of 
anomalous dimensions and therefore allows us to get some insight into the convergence properties 
of the conformal expansion. 

2 Light cone pion distribution amplitudes 

The relation (|2|) of the vacuum-to-pion transition amplitude has to be understood as an operator 
product expansion taken at tree level. At higher orders there will be corrections to eg. the 
coefficient of (j){u), which is taken to unity in (|2|). A more complete form of @ to twist-4 
accuracy is, therefore, (gauge-factors are implicit) 
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where the distribution amplitudes (p{u), gi{u) and g2{u) are defined as matrix elements of non- 
local quark-antiquark operators at strictly light-like separation. The projection on the light- 
cone introduces divergences, and the distribution amplitudes thus have to be renormalised. We 
choose the MS subtraction scheme and the normalisation scale /x^ = 1 /x^ so that the coefficient 
functions are free from logarithms of the type log(x^^'^). 

It is generally believed that the perturbative series r„a" will have factorially divergent 
coefficients This introduces ambiguities in the definition of the sum of the series, usually 

referred to as infrared renormalon ambiguities. These ambiguities are artifacts of the perturba- 
tive factorisation, and for a physically observable transition amplitude these divergences must 
cancel with the corresponding ultraviolet renormalon ambiguities in the definition of the higher 
twist distribution amplitudes. In a different language, the renormalon problem indicates a dou- 
ble counting of the regions of small momenta if perturbation theory is not calculated with an 
explicit IR cutoff, and the double counting is reflected by quadratic divergences of the twist-4 
operators. We refer the reader to Ref. for the thorough discussion. The higher-twist con- 
tributions are, therefore, necessary to render the perturbative treatment of physical observables 
well defined. It follows that the ambiguities of the leading-twist perturbative series can be taken 
as a "minimal model" of higher-twist distribution amplitudes. This we will call the renormalon 
model 

Instead of calculating the higher-order contributions to the leading-twist coefficient func- 
tion directly, it is more instructive to consider the ultraviolet renormalon ambiguity in twist-4 
operators, which is equivalent according to the above argument. The calculation of the ultra- 
violet renormalon ambiguity in the twist-4 operators is similar to a calculation by L. Magnea, 
V.M. Braun and M. Beneke reported in As in the standard approach, it is convenient 
to study first the three-particle quark-antiquark-gluon operators and use the equations of mo- 
tion. Specifically, the following identities, derived from the equations of motion for QCD string 
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operators can be used: 
d 



dx,. 



(i(-x)7^75n(x) = i / dvvd{-x)xagGai3{vx)'yi3jr,u{x), 



df, {d{-x)ji,-f5u{x)} 



i / dvd{-x)XagGai3{vx)'Ji3j5u{x). 



(5) 
(6) 



Here 9^ is the derivative with respect to a translation of the nonlocal operator, while is 
the derivative with respect to the interquark separation. Gai3 is the gluon field strength tensor. 
The three-particle light-cone distribution amplitudes on the right hand side can be parametrised 
as 



di-x)'yf,-f5Ga/3ivx)u{x) 7r+(p) 



1 



px ' " 



-ipx{ai —a2+va3) 



(7) 



+ 



Pl3 9a^^ 



^aPfj, 

px 



Pa g/S^L 



px 



-ipx{ai —a2+va:j) 



The operator on the left hand side contains contributions of both twist-3 and twist-4, but the 
matrix element of the twist-3 part vanishes. The pion decay constant /^r is defined through 



(0|n(0)7^75a!(0)|^+(P)) = i^P^, 
and the integration measure on the right hand side is defined as 



J = J daida2da3(5(l — ai — a2 — as) 



(8) 



(9) 



By sandwiching (^) and (|^) between the vacuum and pion state it is now possible to derive the 
following relations between the parameterisations of the twist-4 distribution amplitudes and the 
three particle operator Q 



52 (n) 
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(10) 



dai / da2— 2 (itai — 'ua2)[2$_L — ^>||](ai, a2, as) 



aq 



Here as = 1 — ai — a2 and u = 1 — u. In the following analysis we will calculate an explicit 
expression for the special combination [2$_|_ — <I>|j](ai, a2, as) in the renormalon model and restore 
gi{u) and g2{u) from these expressions. 



3 The Renormalon Calculation 

It is generally believed that although the perturbative coefficients for QCD perturbative series 
show signs of factorial divergence, the series are asymptotic series for which we can assign a sum 
using the Borel transformation . The Borel transform of a perturbative series 

oo 

R{a) = roa + na'^ + r2a^ + • • • = ^ r„a"+^ (11) 

n=0 
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is defined as 



S[i?](t) = (12) 



n 

n=0 



If B[R]{t) is well-defined on the positive real axis and is integrable we may define the Borel sum 
of the series as 

l-OO 

R= dte'^/^BlRjit), (13) 



with a series expansion identical to that of R. For convergent series, the Borel sum is therefore 
identical to the sum of the original series, and we can regard the Borel sum as an extension of 
the definition of the sum of a series to divergent, asymptotic series. 

However, for a QCD observable R the Borel transform B[R\ has singularities on the real 
axis, leading to ambiguities in the definition of the integral (|l^). These so-called renormalons 
are a direct result of the divergence of the original series, and as already mentioned the residues 
of the renormalon poles are taken as an estimate of the non-perturbative contribution to the 
observable [Q. For QCD observables, one finds that there are infrared renormalons stemming 
from the low- momentum region of higher order loop contributions, situated at w = t/Po = 
1, 2, . . .. Here Pq is the first /3-function coefficient ^^"^^ = —Pocti + 



As seen from (13), the contribution from a pole situated at Wi in the Borel plane to the 
observable will be suppressed by exp{—WiPo/as). As the QCD coupling (at one-loop approx- 
imation) fulfils 1/as = /9o lii((5^/A^), we will need the IR renormalon at if = 1 to obtain the 
l/Q^ twist-4 corrections. 

It is certainly impossible to calculate the residue of the Borel transform at a given pole 
exactly, since that would require full knowledge of all the perturbative coefficients in the series. 
Instead we will apply the so-called Naive non-Abelianization (NNA) |l^, 13] in which the 



coefficients are approximated by their leading /3o part. The NNA starts from the approximation 
of the perturbative coefficients with the part coming from diagrams with one string of fermion 
bubbles inserted to a gluon line. In this way the leading Nj piece {Nf being the number of quark 
flavours) is found. The full one-loop /5 function behaviour is then restored by substituting for 
Nf the scaled /9- function coefficient Nj — ^ . This results in a gauge- invariant approximation of 
the diagrams to all orders in the coupling |^. The positions of the poles in the Borel plane do 
not depend on the higher /3-function coefficients and so are unchanged in the NNA compared 
to the full theory. On the other hand, the residues change and one actually finds that also the 
nature of the divergences may change. In this way, the pole-divergence of the Borel transform 
at tt; = 1 in the NNA is actually the end-point of a branch cut in the full theory. 

The NNA is useful for doing calculations, since in this approximation the Borel transforma- 
tion applies to the modified gluon propagator itself rather than the diagram as a whole. The 
Borel-transformed effective gluon propagator in the Lorentz gauge is given by 

B[a.i,^)D,^]H = ^{9,.-^-^)[-^e~^y , (14) 

where the kf^k^-terms can be neglected for gauge-invariant quantities, fi is the renormalisation 
scale and C is a renormalisation scheme-dependent constant. In the MS-scheme C = —5/3. It 
is noted that the only essential difference between (|l4|) and the original gluon propagator is the 
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Figure 1: The twist-4 contributions to the renormahsation of the three-particle operator (0). 
The third diagram is scale-less and vanishing P]. 



power of k"^ in the denominator, and therefore a calculation of the Borel transform of a diagram 
in the NNA to all orders in the coupling is no more difficult than the calculation of the same 
diagram with a simple gluon propagator. 

We start by noting that to project out the the sought combination [2$_l — $||], we have to 
contract the three-particle operator in (0) with (— x^f/^^). We want to calculate the leading 
ultraviolet renormalon singularity in the matrix element of this operator at w = 1. Since this 
singularity comes from contributions of large momenta, the particular external states do not 
matter and it is sufficient to evaluate the contribution of the diagrams in Fig. || between quark 
states. 

In the Borel representation the matrix element of the projected operator in between 
quark states of momentum p is given by 



P 



d{-xhtil5Gai3{vx)u{x)x°'g''f^ 



-ivkx kp{p + k)a 



21 



where v = 1 — v. For x strictly light-like, we find for the Borel-transformed operators in a 
neighbourhood oi w = 1 

d{-x)j^j5Gap{vx)u{x)x'^g'''^ 



quad, divorg. V ^ ^ ^^ ^2tt'^ W - I 

»1 _ 

daa {d [—x] 759^ u[x{l — av)] — d [—x{l — av)] 75^ n[x]} . (16) 

As already mentioned, the quadratic ultraviolet divergence of the residue of the pole is indepen- 
dent of the external states, which allows us to write the relation as an operator relation. Note 
by comparing ( p!6|) and (|2|) that we have written the twist-4 operator as an integral over leading 
twist contributions. By taking the residue at t;; = 1 of the pion to vacuum transition amplitude 
of the twist~4 operator we find 

[2$^-$l|](a.) = -^^/^'«-''^{(r^^'/'(l-«2)-^^^^<^(ai)}, (17) 
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Choosing the asymptotic expression for the pion distribution amphtude (f>{u) in (|T7|), and insert- 
ing this result in (10) we obtain the following explicit expressions for the distribution amplitudes 
9i{u),g2{u): 



92{u) = N{u\ogu — u\ogu) 
.1 



9i{u) 



N 



N- 



2uu + 2 log u log u + 2u log u + 2u log u 

- 2uL\.2 



— 2 — 2 

-u log u — u log u — 2nLi2 



(18) 



(19) 



2-K 



with Li2(z) the di-logarithm. 



4 Discussion 

Unfortunately, the overall normalisation of the distribution amplitudes gi{u), g2{u) is not acces- 
sible from the renormalon model and has to be found from eg. the QCD sum rule method. This 
is also the case for renormalon models of DIS structure functions and the conformal approach 
to pion distribution amplitudes. 

However, the normalisation factor N found using the renormalon residue can be taken to 
indicate the order of magnitude of the twist-4 correction. Using fi = A = 0.250GeV we get 

~ 0.091GeV^. Thereby we find gi{u)du = 0.0088GeV^ which has to be compared with the 
result obtained by normalising the conformal results using light-cone QCD sum rules. Using the 
functional form and normalisation for the conformal approach of Ref. Q we find gi{u)du = 
0.028GeV^. For the first moment of the distribution amplitude 52 we find du{2u — l)g2{u) = 
O.OOSlGeV^ and /J^dn(2n — l)g2{u) = O.OllGeV^ for the renormalon and conformal approach 
respectively. 

It is more interesting that the shapes of the predicted wave amplitudes are close to the 
ones obtained in the conformal approach of see Fig. |2|. The functions have been normalised 
so that Jq dugi{u) = 1 and du{2u — l)g2{u) = 1. The similarity is encouraging, since it 
indicates that the approximations are well understood. The somewhat broader shape of the 
distribution amplitude gi{u) in the renormalon model compared to the conformal expansion is 
explained by the contributions of higher conformal spin operators. The anomalous dimensions 
of operators are ignored in the NNA approximation, so that it effectively corresponds to the 
summation of the conformal series to all orders, including, however, all kinematic suppression 
factors. The renormalon model, therefore, pretends to make a simple estimate of higher orders 
in the conformal expansion, and the main result of our study is that such contributions seem 
to be comfortably small, apart from in the end-point regions. The asymptotic behaviour of the 
renormalon model prediction for the distribution amplitude (72 (^i) at x — > 0, 1 coincides with 
that of the lowest conformal spin contribution to 52 ('u). At the same time, the behaviour of 
gi{u) is different. The renormalon result behaves as u for u compared to the u'^ behaviour 
of each conformal operator, indicating that the series in conformal operators is not converging 
uniformly. This feature will have an effect on observables sensitive to these regions like i?-meson 
decay form factors 0. 
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gl(u) g2(u) 




Figure 2: gi{u) (left) and g2{u) (right) calculated in both the renormalon model (line) and 
the conformal wave expansion of |^ (dots for the lowest conformal spin result and dashes for 
the result including next-to- lowest conformal spin operators) . The distribution amplitudes have 
been normalised so that dugi{u) = 1 and du{2u — 1)5(2 (n) = 1. 
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